function y=Xa(x)
    y=exp(-3.837)*x^0.602;
end